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SOME  ORBITAL  ELEMENTS  USEFUL  IN  SPACE 
TRAJECTORY  CALCULATIONS^ 

William  Kizner 

Jet  Propulsion  Laooratory 
California  Institute  of  Technology 
Pasadena,  California 

ABSTRACT 

A  set  of  orbital  elements  is  described  which  is 
applicable  when  the  motion  is  hyperbolic  and  even 
rectilinear,  such  as  the  motion  associated  with  a  lunar 
landing.  This  set  provides  a  convenient  description  of 
miss  distance  and  furnishes  simple  relationships  for 
many  complex  geometrical  problems  in  lunar  and  inter- 
planetarv'  flight.  The  use  of  these  elements  with  others 
which  are  described  enables  the  computation  of  a 
variety  of  differential  corrections. 


^This  paper  presents  the  results  of  one  phase  of  research  carried  out  at 
the  Jet  Propulsion  Laboratory,  California  Institute  of  Technology,  under  Contract 
No.  NASw-6,  sponsored  by  the  National  Aeronautics  and  Space  Administration. 
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I.  INTRODUCTION 

Space -flight  research  characteristically  involves  the  calculation  of  innumer¬ 
able  trajectories,  some  of  which  are  reference  trajectories  and  others  variations  of 
reference  trajectories.  The  orbital  eleme.nts  given  here  can  be  used  to  determine 
differential  corrections  or  changes  from  the  standard  trajectory  without  the 
necessity  of  coinputing  varied  trajectories,  providing  a  welcome  reduction  in 
computation  time. 

As  an  e.xample  of  their  use,  these  orbital  elements  can  be  applied  in  the 
determination  of  a  trajectory  which  impacts  the  target  in  a  given  way.  If  the  target 
were  on  the  Earth's  surface,  the  miss  distance  would  be  described  in  terms  of 
range  and  azimuth  error,  and  a  typical  differential  correction  would  be  the  partial 
derivative  of  range  error  with  respect  to  speed  at  burnout.  However,  if  the  target  is 
the  Moon  and  a  vertical  impact  is  desired,  some  other  measure  of  miss  distance  has 
to  be  used.  By  analogy  with  the  terrestrial  target  calculations,  two  components  are 
necessary.  For  this  purpose,  with  a  lunar  target,  use  is  made  of  two  components  of 
B,  the  impact  parameter  vector,  which  has  magnitude  b  and  lies  in  the  orbital  plane. 
It  is  directed  from  the  focus  perpendicularly  to  the  incoming  asymptote  of  the 
osculating  hyperbola,  evaluated  at  a  time  when  the  vehicle  is  near  the  target.  The 
elliptic  case  is  not  treated. 

Having  described  the  miss  distance  in  terms  of  orbital  elements,  one  can 
obtain  analytically  the  approximate  changes  in  these  elements  due  to  changes  in  the 
initial  osculating  conic  referenced  to  the  Earth.  Of  course,  it  is  not  necessary  to 
obtain  these  changes  analytically  in  order  to  use  the  elements  as  a  convenient 
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moasuro  of  miss  dislancf.  When  it  is  known  how  tin-  miss  emeponems  vary  with 
changes  in  the  initial  conditioiis,  either  approximately  or  by  a  numc'rical  solution 
(jf  the  varic'd  equations  of  motion,  it  is  jtossible  to  arrive  at  anv  desiri  d  kind  of 
impact. 

In  this  paper  the  general  theory  and  metliods  of  computation  are  given. 

The  results  will  be  diseussiui  in  subsequent  pape-rs.  The  use  of  orbital  elements 
for  diffeia.mtial  corrections  is  described  by  Hckert  and  Hrouwer  (Ref.  1)  and 
S.  Merrick  with  otheu'S  of  Aeronutronic  (Ref.  2). 
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II.  THE  B  SET  OF  ELEMENTS 


For  near -rectilinear  hyperbolic  motion  the  B  set  of  elements  has  been  found 
to  be  a  well-behaved  set  of  parameters  that  produces  no  singularities,  provided  th(? 
motion  does  ?int  become  pai'abolic.  The  variations  of  this  set  with  respect  to  changes 
in  the  initial  coordinates  are  surprisingly  linear.  This  set  is  being  used  e.'iclusively 
at  the  Laboratory  to  desigti  trajectories  (Ref.  3). 

The  B  vector,  whic !i  fornis  tlie  basis  cif  the  B  set,  is  defined  as  a  vector 
originating  at  the  focus  of  hyperbola  and  directed  to  the  incoming  asymptote  (see 
Fig.  1).  Analytically, 

B  =  b^J^P  .  f:Q'j  (1) 

where 


with  cj  the  angular  moinentuni  <  onstant,  .  •  ihe 
and 


(21 


is  vuva.  or  twice  tlie  total  energy; 


(3) 


The  other  symbols  have  the  usual  meaning.  When  defined  in  this  way,  B  can  be 
thought  of  a.s  the  vector  miss  which  would  occur  if  the  target  had  no  mass,  although 
this  is  not  strictiy  correct  bec;ause  the  mass  ol  the  target  does  influence  the  o.sculating 
elements. 
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The  dii-cction  of  the  inconiiriff  asymptote  is  specifieti  bv  the  unit  vector  S: 


S  =  -  ?  7  t  '^Q 

(•  c 


(4) 


It  can  be  shown  tliat  in  attempts  to  hit  a  farf;et,S  will  liave  approximately 
the  same  orientation  for  all  search  ti-a jectoi-ies.  Since  B  is  perpendicular  to  S, 
it  follows  that  B  can  be  I'epresemed  approximately  Ity  two  numbers,  or  compo¬ 
nents.  The  components  chosen  are  aion^^  unit  vectors  T  and  R,  perpendicular  to  S 
and  to  each  other.  T  lies  in  the  eouatoriai.  ecliptic,  or  other  convenient  fixed 
reference  plane,  and  indicates  uhether  the  miss  is  to  the  right  or  left  (direct  or 
retrograde).  'Flie  only  condition  for  the  fixed  reference  plane  is  that  it  is  not 
perpe-ndicular  to  the  asymptote.  The  other  vector,  R,  indicates  the  up  or  down 
component  of  miss.  The  two  cf)mponents  are  H  •  T  and  B  •  R.  called  m^  and  m2- 
If  S  is  spt.-<.'ifi(.'d  bv  ;wo  angl>.-s  <'ind  defined  bs’ 


sin  ^ 


os  cp^ 


os  I  OS  cp^  =  S._ 


sin  cos  <p^  = 


(5) 

(G) 

(7) 

(8) 


tlien  T  is  giveti  bv 


T^  =  sin  Og 


T^,  =  -  cos  dg 


(H) 

(10) 
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(11) 
(12) 

A  positive  value  for  B  •  T  indicates  that  the  motion  is  direct;  a  positive  value  for 
B  ■  R  indicates  that  the  motion  is  above  the  target. 

If  the  approximation  is  used  that  S  is  constant  for  a  given  class  of  trajec¬ 
tories,  and  that  the  components  m2  and  m.2  are  linear  functions  of  initial  errors, 
tlien  valuable  information  can  be  obtained  analyttcally  bv  calculating  onlv  one 
trajectory  near  tlie  iai'gei.  For  exam[)le,  it  is  possible  by  a  relatively  simple 
proc  edure  to  find  a  mapping  between  the  miss  components  and  the  selenographic 
coordinates.  Since  ’he  probability  distribution  for  impacting  can  be  easily  deter¬ 
mined  in  terms  of  the  rtiiss  components  (where  linear  theory  is  applicable),  the 
distribution  in  terms  of  the  selenographic  coordinates  ntay  be  obtained  by  the 
known  mapping. 

The  mapping  miss  compon'‘“nts  to  selenographic  coordinates  depends 

e 

essentially  on  two  relations,  shown  in  Fig.  2,  The  line  OP  is  the  path  directed  at 
the  center  of  the  target.  Tliis  line  also  represents  the  asymptote  of  the  osculating 
conic,  which  is  very  nearly  parallel  to  the  asymptotes  of  other  conics  which 

9 

intersect  the  surface  of  the  target.  Consideration  is  now  given  to  the  case  where 
the  impa' t  is  not  vertical  and  the  impact  parameter  is  along  the  T  axis.  Since 
the  plane  ol  the  trajectory  contains  the  original  S  vector,  the  loci  of  the  impact 
points  on  the  surfai  e  of  the  target  (con.sidered  a  sphere)  must  he  a  great  circle. 
This  is  true  whenever  the  direction  of  B  is  lield  fixed  and  ‘.he  magnitude  is  varied. 
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Tiiis  i.s  the  fii's',  relation.  The  othoj-  is  that  the  anf^le  between  ^reat  cii’cles  at  P 
is  the  same  as  the  ancle  bet  weim  the  eori'espondinf^  B  vectors. 

Ancjther  oxtimple  of  tlie  application  of  this  method  is  in  searching  for  a 
ti  ajectoi  N  whi(  h  p.asses  ov’cr  the  pole  of  the  Moon.  The  conditions  can  be  con- 
v'ertc-d  into  miss  l  omponenis  by  fii'st  speciuing  the  pericenter  distance  q;  this  is 
equi'/alcm  to  specifying  the  altitude  at  closest  approach.  Then  the  required  value 
of  b  is 

b  =  -vy/q“  -  2aq  (13) 

where  a  is  negative  for  tlie  hyperbola. 

Since 

1  /  2  2 

b  -  ,/  nil  ^  ”'2  (14) 

tin;  ratio  o:  the  components  and  their  signs  must  In;  determined.  This  is  found  b>’ 
projee.ing  a  unit  v('('’or  in  the  direetion  of  tlie  pole,  O,  on  the  T  and  R  axes.  Then 

_ _ 

b  y(o  •  T)“  +  (O  •  K)“ 

^^^2  _ _ _  (16) 

b  yio  *  f)^  +  (O  •  R)^ 

Thus  a  seemingly  complex  problem  is  solved  by  the  use  of  three  simple  equations, 
(13),  (15),  and  (IG). 
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Tht,'  .six  orbital  elements  that  are  used  are  the  epoi  h  of  periapsis  T,  mj,  m2, 
‘^3’^s’  '^S‘  find  the  rectangular  coordinates  from  these  elements,  Cj  is 

found  from  (14)  and  (2),  S  from  (5),  (7),  and  (8),  T  from  (9)  to  (11),  and  R  from  (12), 
from  which 

B  -  ni]  T  +  m2  R  (17) 

From  (1)  and  (4),  P  and  Q  are  found.  If  the  epoch  is  j^iven,  Kepler's  equation  can  be 
solved  to  find  F,  and  the  rectangular  coordinates  may  be  computed.  Special 
techniques  foi  s.Mviti^  Kepler's  e<]uation  have  been  developed  v/hich  apply  when  the 
eccentricity  is  eiose  to  1. 
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► 


III.  THE  w  SET  OF  ELEMENTS 


When  the  orbital  plane  is  well  determined,  the  w  set  of  elements  is  used:  T, 
Cj,  Cg,  w^,  W2,  and  Wg.  The  angles  represent  rotations  about  the  P.  Q,  and  W 
axes  in  a  given  order  and  have  been  defined  so  that  they  apply  to  finite  rotations. 

Here  it  will  be  assumed  that  the  rotations  are  infinitesinual. 

The  use  of  the  rotational  elements  is  similar  to  that  described  by  Eckert  and 
Brouwer  (Ref.  2),  who  use  rotations  about  the  P,  Q,  and  W  axes.  The  rotations 
Wp  W2.  and  W2  are  about  the  reference  P,  Q,  and  W  axes.  The  other  elements  were 
chosen  so  that  they  would  apply  to  elliptic  motion  with  a  large  eccentricity,  to 
parabolic,  and  to  hyperbolic  motion.  For  small  eccentricities  where  the  periapsis 
is  not  well  determined,  -nT  +  w^  may  be  used  as  a  variable,  similar  to  Eckert  and 
Brouwer.  Since  these  elements  must  apply  to  motion  where  the  eccentricity  is  close 
to  1,  it  was  decided  to  eliminate  the  eccentricity  as  a  variable  because  it  is  not  well' 
behaved.  Also,  a  catmot  be  used  because  it  is  discontinuous.  Instead, 


'^3 


h 

a 


(18) 


is  used,  and  is  well  behaved.  The  reason  for  the  choice  of  T  instead  of  M  as  an 
element  follows  if  the  parabolic  case  must  be  inc  .’uded.  Thus  T,  c^,  c^,  and  the  wo 
w'ere  chosen  as  orbital  elements. 

These  elements  are  used  in  two  different  ways  depending  on  wdiether  the 
are  finite  (as  in  guidance  studies)  or  infinitesimal  (as  in  differential  correciion 
computations).  If  the  w^  are  to  be  kept  small,  the  standai'd  P,  Q,  and  W  are 
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taken  to  be  the  osculating  P,  Q,  and  W  along  the  reference  tiajectory. 
only  are  the  variations  small,  but  the  rotations  are  about  "body- fixed" 
Li  other  words,  if  initially  a  rotation  is  performed,  at  a  latei-  time 
represent  a  rotation  about  the  P  axis  of  the  later  osculating  conic. 


Then  not 
a.xcs. 
this  will 
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IV.  DIFFERENTIAL  CORRECTIONS 


A.  Orbital-to -Rectangular  Changes 

The  formulas  whic'h  give  the  changes  in  inei'tial  rectangular  cooi'dinates 
due  to  changes  in  the  orbital  elements  will  now  be  presented. 

1.  Chatiges  in  the  w  set.  The  w  set  is  treated  first,  since  some  of  its 
relations  are  further  utilized  in  the  transformation  of  the  H  set. 

The  equations  for  the  hyperbola  read 

r  =  aP(cos!i  F  -  e)  +  bQ  sinJt  F  (L'O 

V"  ! —  —  sinh  F  ,  , —  Qcosh  F  /on^ 

r  =  P  — -  t  ypp  - -  (20) 

Differimtiating ,  ;ind  kccpiiie  l.  ihe  cp-''  !!  ai  whit-h  the  coordinate.'S  arc  ob.servcd, 
fixed,  \'ield^ 


dr  =  daP(cosh  F  -  e)  *  dPa(cosh  F  -  e) 


+  dFaP  sinh  F  -  deaP  r  dbQ  sinh  F 


(2  1) 


+  dQb  sinh  F  +  dFbQ  cosh  F 


dr  =  da  ,7 


q  sinli  F  —  /  sinh  F 

P - p -  -  dP  s/-Ata  P  - - 


ai-  / -  D  co.sh  F  ,  ~  sinh  F 

-  dl-  V -/ua  P  — —  —  +  dr  v/ -qa  P  - - — 

r  1-^ 

+  dp  A  yj  +  dqy^t 


(22) 


+  dF  y  up  Q  --(AL-  -  dr  y/Ttp  Q  — 


cosh  F 
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These  equations  as  they  stand  are  not  in  terms  of  the  orbital  elements  T, 
Cl,  C3,  wi,  W2,  and  W3.  Each  will  be  expajided  in  tui'n.  From  the  definition 


da  =  - 


—  dco 
Co  a 


(2  3) 


dP,  dQ,  and  dW  ai'e  found  from 


dP  =  ”  Ww'2 

dQ  =  Ww  |  -  Pw'3 

dW  =  Pw2  ■  Q^'  1 


(24) 

(25) 

(26) 


as  can  be  seen  from  drawinji  a  simple  picture  of  any  rotation.  Equations  (24)  to  (26) 

o 

can  be  put  in  matrix  form: 

hx  l^x  «x  'Vx\  /  »  -»3  ''■2\ 


dP.^.  dt^^.  dW^. 


w  3  0  -  w 


(27) 


P,  0,. 

dW,y  \p,  W.,y  wi  0  I 

This  form  is  suitable  if  the  rotations  are  finite,  and  the  last  matrix  represents 
the  derivatives  of  an  orthogonal  matrix  representing  the  rotation.  Also  by  direct 
differentiation, 


di; 

dh 

dn 

dp 


2cp2  (  1 

dc  1 


1  JL 

2  03 
2  c.  jdc- 


dc 

1 


3 


(28) 

(39) 

(30) 

(31) 
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Inc  idenialK  ,  (23)  lu  (31)  also  apply  when  tlie  conic  is  an  cTiipse,  Keplei’'s 
equation  is  differentiaU'd: 


ciK 


dn(  t  -  T)  -  ndT  -  de  sinh  h' 
c?  I  fcsn  F  -  1 


(32) 


where  dn  and  de  have  alreadv  been  j^dvon  in  terms  of  the  oi'bital  elements.  Also, 


dr  =  -d3(t-  cmsli  h'  -  1)  -  clc-  a  oosh  F  -  dFae  sinh  F  (33) 

The  formulas  for  the  elliptic:  case  arc-  so  similar  that  they  are  not  listed. 

2 ■  Chan^u-s  in  the  B  set.  For  this  set,  liie  division  between  parameters 
which  delei'mine  the  orientation  (the  wp  and  parameters  which  detei’mine  the  other 
propej'ties  is  lost.  Mc;st  <>■  tlie  elements  affcc  t  both  the  orientation  and  shape  of 
tlie  conic.  However,  Fiqs.  (21)  and  (22)  hold.  Since 


P 


(34) 


and 

o 


(33) 


the  variations  in  B  and  Q  c:an  be  dc-rived  if  the  right-hand  side  is  differentiable.  For 
the  rectilinear  case,  when  !c  a;jproac:hes  zero,  the  term  in  (35)  involving  B  can  cause 
difficulty.  Since  the  rectilinear  ca.se  presents  .some  analytical  problems,  it  ■v\ill 
bo  fully  treated.  The  nonrectilinear  case  iS  straighiforwai'd.  and  will  not  be 
analyzed  in  detail. 
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From  Eqs.  (28)  and  (31).  de  and  dp  are  /,ero  siiu  <  i  ^  /<  ro  :  jc  th*  r  ec  t: 

linear  case.  Hence,  the  equations  of  motion  (19)  and  (20)  can  be  differentiatt'd 
keeping  e  =  1  and  p  =  0: 


dr  =  daP(cosh  F  -  1)  +  dPa(cosh  F  -  1) 


+  dPaP  sinh  F  f  dbQ  sinh  F 


(36) 


1  /  q  sinh  F 


dr  =  da  -  /-  T  P 

■>  V  a 


-  dP 


pa 


sinh  F 


(3  7) 


OF  y^a  P 


—  ^  cosh  !•' 


+  dr  -y-pa  P 


1- 


sinli  F 

.2 


It  should  be  notec  that  Q  appears  in  Eq.  (36)  and  has  not  been  defined  for 
the  rectilinear  case.  It  will  be  d(.>fined  so  that  it  depends  on  each  variation  in  a 
natural  way.  First  da  and  OF  are  given  by  E()S.  (23)  and  (32)  as  before.  Differ¬ 
entiating  (34)  and  evaluating  it  for  the  rectilinear  case, 

dP  =  dS  r  dB  ^  (38) 

d(^'  i s  not  used . 

dS^  =  -  'kin  dg  cos  <^g  ddg  -  cos  ''g  sin  c^gdc^:^ 

dSy  =  fos  dg  cosc^gdOg  -  sin  Og  sin<f>g  d<pg  (39) 

dS^  =  cos  cpg  d<^g 

dB  =  d(I3  •  T)  T  +  d(B  '  R)  R  (40) 


T  and  R  arc  given  una mlhguouslv  by  F<is.  (a)  to  (12). 
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It  now  reir.ains  to  define  Q,  Since  dB  is  contained  in  tlie  T,  R  plane  from  (40), 
it  IS  perpendicular  to  S,  or  to  F  which  is  equal  to  S.  Hence,  there  is  no  reason 
why  Q  should  not  be  defined  to  he  in  the  opposite  direction  to  dB.  This  will  insure 
that  B  is  a  continuous  function  for  the  variation  in  any  element.  Hence, 


Q  = - 


m2  +  ni2 
1  2 


'^m2  +  m2 

1  2 


Since 


'm2  +  m2 
1  2 


m  j  m^ 

IT  "^^2 


This  formula  must  be  modified  when  b  goes  to  zero,  which  is  the  rectilinear  case. 
If  it  is  as.sumi'd  that  dm^  or  drn2  is  positive,  then 


db  =  dm  2 


for  an  variation,  and 


db  =  dmr 


for  an  m2  variation.  In  general. 


db  =\/{dmi)^  +  (dm2)^ 


!'nr  the  rectilinear  case,  c^  is  no  longer  an  independent  variable 


dc^  =  db 
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B.  Rectangulai  -lo-Orbilal  Changes 

1.  Changes  in  the  w  set.  It  is  a  simple  matter  to  find  the  changes  in  the 
angular  momentum  and  energy  constants  due  to  changes  in  the  rectangular 
coordinates: 


(48) 

(49) 

(50) 


For  the  rectilinear  case  W  is  detei  niined  by  the  perturbation,  and  Eq.  (44) 
reduces  to 

dc][  =  Idr  X  r  +  r  X  dr  1  (51) 

which  can  also  be  derived  from  the  fact  tha*  Cj  originally  is  zero.  Since 


.  .  2/j 

Co  =  r  *  r  -  — I 

1  j’  I 


dc3 


(52) 

(53) 


The  formulas  for  changes  in  a,  e,  p,  and  n,  which  depend  only  on  c^  and  cg,  are 
unchanged.  For  the  hyperbola  from  (33), 


dF  =■ 


+  da(e  cosh  F  -  1)  +  de(a  cosh  F) 
ae  sinh  P' 


(5  j) 
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(IT  =  -  ^  de  sinli  F  +  dF(e  cosh  F  - 

1)  -  dn(t  -  T)j  . 

(56) 

Tlie  changes  in  orientation  will  next  be  given. 

W  = 

^'1 

Since 

(57) 

dW  =  w 

(  1  Cl 

dci 

•  (58) 

From  Eq.  (2  6), 

dW  •  P  =  W2 

(59) 

d\V  •  Q  =  -wi 

(60) 

A 

To  find  Wg,  v  is  introduced: 

A 

V  =  V  +  W3 

(61) 

A 

The  term  v  represents  the  polar  angle  that  the  vehicle  makes  with  the  reference 
P  axis  for  infinitesimal  l  otations  (see  Pig.  3). 

Hence, 

A 

dw3  =  dv  -  dv  (62) 
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Since 


cos  V 


\ 


(63) 


dv 


_ 1__ 

sin  V 


de  !  p 

^Vr 


er 


+ 


dr  p 
er^ 


(64) 


A  —  P  sin  V  Q  cos  V 

dv  =  dr  •  -  -I-  - 

L  -r  r  . 


(65) 


This  last  formula  can  easily  be  derived  by  noting  that  dr  •  P  and  dr  '  Q 
are  the  components  of  dr  on  the  P  and  Q  axes.  A  change  in  position  along  the  P 
axis  results  in  a  rotation  of  -  sin  v/r,  etc. 

Thus  the  changes  in  the  w  set  can  be  computed. 

2.  Changes  in  the  13  set,  rectilinear  case.  The  changes  in  c  C;^,  T,  etc. 
remain  the  same  as  in  the  w  set,  but  the  changes  in  nn,  m2.  <ps<  Og  must 
be  computed.  It  is  convenient  to  define  a  unit  normal  to  W,  for  each  variation  in 
the  rectangular  coordinates.  Since 

cjW  =  r  X  r  (66) 


and  dW  is  defined  to  be  zero. 


dcjW 


dr  X  r  -t  r  X  dr 


(67) 
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or 

dr  X  !•  1-  r  X  dj- 

W  =  ,  -I 

(dr  X  r  t  ;■  X  dr  ( 

(o:)) 

As  an  example  for  a  variation  in  x, 

U-  =  r 

(fid) 

/to  -  O  :  / 

7  ^  .  j 

\  ^  / 

B  -  b(S  X  W) 

(70) 

dH  =  dlj(S  X  W) 

(71) 

d(B  •  T)  -  (dB)  •  T 

(7.') 

d(H  •  K)  --  (dB)  •  K 

(73) 

Foi'  the  rectilinear  or  alinost  l  ectihuoar  <  ase, 

“*■  » • 

S  ^  T 

s 

(74) 

where 

(7.7) 

This  assumes  that  the  vehicle  is  a  .sufficient  distance  fcom  the  focus  that  it  i.s 
moviiif^  parallel  to  the  asymptote.  In  p-  attice,  tlrj.  (74)  will  be  an  (excellent 
approximation  if  the  otlirr  conditions  for  the  linear  c;jse  aie  satisfied.  Hence. 
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dr  T  V  •  di' 

db  =  ~  -  b  — “n  — 


(70) 


From  (39), 


d(/)g 


dSz 
’  OS(^S 


dSy  T  sin  Og  sinc^g  dt^g 
eos  dg  coS(^g 


(77) 

(78) 


C.  Differentia!  Correotion.s 

In  the  general  problem  (jf  differentia !  c'or rectic^n.s  one  is  eonr  erned  with  the 
changes  in  the  position  and  velocity  eo(mdinates  of  sonie  point  along  the  path  due  to 
changes  in  the  orbital  elements  at  a  different  point.  For  the  two-body  problem  the 
path  can  be  represented  analytically  and  the  deviatioits  can  be  obtained  rigorously. 
Even  if  the  path  is  n<d  e.xactly  an  ellipse,  it  may  be  possible  to  eaiculate  changes  in 
the  path  as  if  a  conic:  solution  applied.  The  free  flight  paili  of  an  ICBM  diffc'rs  so 
little  from  an  ellipse  that  in  this  case  the  deviations  may  be  computed  on  tlie  basis  of 
formulas  that  are  derived  from  the  two-body  probUnn. 

This  program  of  m'ror  correction  can  be  used  only  in  tiie  cast*  whei’c  the 
actual  path  is  almost  a  conic  section  and  the  elements  chosen  exhibit  the  ftnlowing 
propertie.s.  Consider  a  reference  trajectory  at  tiir.e  t^.  One  of  its  elements  vp-  (ty) 
is  subjected  to  a  small  perturbation  (tg),  while  all  the  other  elements  are  kept 
fixed.  One  now  computes  the  path  from  tg  to  a  later  time  tj.  denoting  the  elements 
of  the  varied  path  by  q,;.  If  ik^w 

di  (ti)  -  di  (tl)  ==  i-li  (tg)  i  =  I--  -  6  (79) 


1 
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► 


and 

"qj  (tj)  -  Mj  (ti)  -  0  i  j  (80) 

then  the  elements  are  well  behaved  and  no  special  variational  equations  need  to  be 
solved  since 


<5qi  (to) 


(81) 


The  orbital  elements  that  have  been  chosen  exhibit  this  property  for  trajectories 
encountered  in  deep  space  exploration  where  eccentricities  are  not  close  to  zero.  It 
has  been  shown,  however,  that  a  lunar  or  interplanetary  trajectory  can  be  broken  up 
into  segments  each  of  which  is  in  a  region  of  a  dominating  body.  Let  us  examine  in 
this  light  a  trajectory  which  is  meant  to  intersect  with  the  Moon  and  calculate  the 
differential  corrections  of  errors  at  the  Moon  as  a  function  of  errors  at  injection. 

Using  the  w  set  at  injection  (the  orbital  plane  being  well  defined),  one  converts 
the  various  errors  at  injection  into  orbital  elements  using  the  osculating  conic  sections 
at  injection.  Since  the  region  of  influence  of  the  Moon  is  a  concentric  sphere  of 
66,000  km  radius,  this  is  taken  as  the  boundary  between  Earth  and  Moon 
domination.  At  this  point  the  changes  in  the  elements  of  the  osculating  conic  sections 
referenced  to  the  Earth  have  to  be  changed  to  those  referenced  to  the  Moon.  This 
transformation  is  accomplished  in  three  steps.  The  first  step  is  to  convert  the 
orbital  changes  of  the  Earth -referenced  osculating  conic  into  variations  in  an  Earth- 
referenced,  nonrotating,  cart<!sian  coordinate  system.  The  next  step  is  to  convert 
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the  cartesian  coordinate  changes  to  a  Tvloon -j  eferenced  cartesian  system.  If  tiie  two 
coordinate  systems  are  parallel,  it  follows  tiiat 

^12  +"^2  "  '’1 
kl2  +  i'2  "  "’I 

where  R]^2  position  of  the  second  center  cm"  attraction  with  respect  to  the  first, 

r2  is  the  position  of  the  vehicle  with  respect  to  the  second  center,  r^  is  the  position 
of  the  veliicle  with  respect  to  the  fii'St  center';  -nui  similar!',  for  the  velocities.  Tims 
the  variations 

or^  =  hr  2 

are  continuous  at  the  boundary  between  the  two  regions.  The  third  step  is  to  convert 
the  Moon-i-eiitered  cartesian  system  int>-  clianges  in  the  osculating  conic  referenceu  t  - 
tile  Moon.  Tlie  osculating  conics  used  in  the  calculations  are,  of  course,  the  ones  at 
the  boundary  and  consequently  may  differ  materially  from  the  initial  as  weii  as  from 
the  final  values. 

In  the  neighborhood  of  the  tai'get  the  B  set  is  used.  With  this  analytical 
method  one  obtains  the  partial  derivatives  of  the  miss  roinponeriUj.5^t  the  taiil^iet  with  a 
precision  of  a  few'  per  cent.  7’he  errors  incuri-ed  are  due  mainly  to  effects  a:  the 
boundary  of  the  tw'o  regions.  Studies  are  under  way  to  circumvent  this  problem. 
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V.  CONCLUSIONS 

The  B  set  and  the  w  set  have  proven  to  be  very  useful  in  the  determination  of 
differential  corrections.  The  B  set  can  be  used  for  cases  where  the  motion  is 
hyperbolic, and  the  w  set  is  used  for  all  motions  except  rectilinear  motion  and  nearly 
circular  motion.  The  combination  of  the  two  sets  may  be  employed  to  predict  changes 
in  the  final  conditions  due  to  changes  in  the  initial  conditions.  Also,  corrections  in 
the  orbit  can  be  made  by  the  same  technique. 

The  method  of  solving  the  ciuations  of  motion  by  variation  of  parameters  may 
be  derived  from  the  ec}uations  given  here.  It  is  possible  to  obtain  the  t'ffects  of  the 
perturbations  on  the  two-body  problem  by  treating  the  perturbations  as  incremental 
velocities.  From  the  transformation  of  velocities  ro  orbital  elements  (presented  here 
for  certain  cases)  the  equations  for  the  variation  of  parameters  follow. 

Anotlier  application  of  these  elements  is  in  guidance  studies.  Here  the  w  set 
is  employed  for  finite  rotations.  For  this  application  it  is  usefui  o  have  the  end 
conditions  in  terms  of  the  B  set  represented  by  a  polynomial  as  a  function  of  tlie 
initial  orbital  elements  of  the  w  set.  Here  the  varied  traiectories  are  calculated 
numerically  to  find  the  polynomial. 

The  complete  equations  used,  including  the  method  of  solution  of  Kepler's 
equations,  will  be  given  in  a  forthcoming  report. 
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a 


A 

b 

B 

c 

'■l 

"3 


e 


E 


F' 


nip 

n 

O 

P 

P 

q 

Q 

I' 

7 

R,  T 


NOMENCLATURE 

semi  major  axis 

rotation  matrix  defined  by  Eq.  (20) 
semiminor  axis 

the  impaet  parameter  vector  defined  by  Eq.  (1) 
defined  for  the  hyperbola  by  Eq.  (3) 

2 

the  angular  momentum  constant  =  r  v 

2 

the  \'is  viva,  or  total  energy  =  s  -  (2p/r) 
eccentricity 

eccentric-  anomaly  for  an  ellipse 
corresponding  anomaly  for  the  hyperbola 
miss  components  =  B  •  T  and  B  •  T,  respectively 
mea.u  motion 

unit  vector  in  the  direction  of  the  Moon's  pole 

semilatus  rectum  of  the  conii. 

unit  vector  directed  toward  the  periapsis 

distance  from  focus  to  vertex  of  conic- 

unit  vector  in  orbital  plane  normal  to  P 

distance  from  focus 

position  vector 

unit  vectors  perpendicular  to  S  which  are  associated  with  the  miss 
components,  defined  by  Ecjs.  (9)  to  (12) 
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NOMENCLATURE  (Cont’d) 


s 

S 


t 

T 

V 

A 

V 

Wi,  W2,  W3 

W 


<^g,  (1^ 


the  speed 

unit  v'ector  in  the  direction  of  the  incoming  asymptote. 

defined  by  Eq.  (4) 

time,  epoch 

epoch  of  periapsis 

true  anomaly 

<it‘iine<i  by  E().  (01) 

establish  the  rotation  of  the  P,  Q,  and  W  axes,  defined  by  Eq.  (20) 
the  unit  vector  ncmmal  to  the  plane;  used  with  P  and  Q 
gravitational  lonstan'  multiplied  by  the  mass  of  the 
attracting  centei' 

angles  which  give  orientation  of  S,  defined  by  Eqs.  (5)  to  (8) 


25 


JPL  Technical  Release  No.  34-84 


Fig.  1.  Geometry  of  hyperbolic 
path  near  target 


Fig.  2.  Effect  of  variation  of  B  •  T  and  B  •  R 
on  surface  of  target 


VARIED  POSITION 


Fig.  3.  The  relationship  of  v,  v. 
and  Wg 


REFERENCE  POSITION 

t  'V, 
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26 


i:lki-:re.\(.'ks 


Eckei’t,  W  .  J.  and  Brouwoi',  D  ,  "The  I  se-  ui  HuctLi 
the  Difier'ontial  Corr'ection  oi  O.'bits,  '  Ast  i  (>nomical  .Trj’v 
No.  13  {.August.  Hi,  1937).  p.  125. 


Coordinates  in 
irr-.al  ,  \'(>L  IvL",  i 


•Aeronnt ronic  Systems,  Inc.  .  Space  \h.-hif.lo  ilphemeris  and  Differential\ 

(  on  (-ction  Program  --  Unified  Thi'orv,  .-\e romitronic  Publication  U-908;' 
Newport  Roach,  California,  Jun<-  14,~l9(i() 

A 

Kiznei  ,  W  .  ,  A  Alethod  of  Describing  Aliss  I)istuir-es  tor  Lunar  and  Inter- 
r J'ajecturies.  External  Publication  No.  674,  Jet  Propulsion 
Laboiatot’/,  lasadiina,  C'alifornia.  Augu.st  1,  li,'59. 


